home *** CD-ROM | disk | FTP | other *** search
/ NetNews Offline 2 / NetNews Offline Volume 2.iso / news / comp / lang / c-part2 / 10668 < prev    next >
Encoding:
Internet Message Format  |  1996-08-05  |  1.2 KB
  1. Path: news2.noc.netcom.net!news
  2. Newsgroups: comp.lang.c
  3. Subject: Re: goto
  4. Message-ID: <314E0BDE.7933@willows.com>
  5. From: Tarang Deshpande <tarang@willows.com>
  6. Date: Mon, 18 Mar 1996 17:20:30 -0800
  7. References: <Pine.OSF.3.91.960313102715.10701D-100000@io.UWinnipeg.ca>
  8. Organization: NETCOM Network Operations
  9. NNTP-Posting-Host: daffy.willows.com
  10. Mime-Version: 1.0
  11. Content-Type: text/plain; charset=us-ascii
  12. Content-Transfer-Encoding: 7bit
  13.  
  14. Bill Simpson wrote:
  15. > 
  16. > There was a goto thread lately, and my problem is to state this algorithm
  17. > cleanly without gotos (which I think is easy, but my attempts have been
  18. > failures).
  19. > 
  20. > 0. x=0;
  21. > 1. x+=erand(maxmean);
  22. > 2. if (urand2()>rate(x)/maxrate)
  23. >         goto step 1
  24. > 3. if (x<=XMAX)
  25. >         {
  26. >         setdot(x,y,z);
  27. >         goto step 1
  28. >         }
  29. > explanation:
  30. > erand(mean) returns exponential double random number with given mean
  31. > urand2() returns uniform double random number between 0 & 1
  32. > rate(). well the rate of the Poisson process varies with x
  33.  
  34.  
  35. x = 0;
  36. do
  37. {
  38.     double uRand;
  39.     double rateCalc;
  40.  
  41.     x += erand ( maxmean );
  42.     uRand = urand2 ();
  43.     rateCalc = rate ( x ) / maxrate;
  44.     if ( urand <= rateCalc )
  45.         if ( x <= XMAX )
  46.             setdot ( x, y, z );
  47. } while ( ( x > XMAX ) || ( uRand > rateCalc ) );
  48.  
  49.